Scaling of the anomalous Hall effect in SrRuOs 
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We measure the anomalous Hall effect (AHE) resistivity p^y in thin films of the itinerant fer- 
romagnet SrRuOs. At low temperatures, the AHE coefficient Ra varies with p^^; ,and at higher 
temperatures, Rs reaches a peak and then changes sign just below Tc- We find that for all films 
studied Ra scales with resistivity in the entire ferromagnetic phase. We attribute the observed be- 
havior to the contribution of the extrinsic side jumps mechanism and the intrinsic Karplus-Luttinger 
(Berry phase) mechanism including the effect of finite scattering rates. 

PACS numbers: 75.47.-m, 72.25.Ba, 75.50.Cc, 72.15.Gd 



I. INTRODUCTION 

Being one of the most intriguing manifestations of a 
transport phenomenon that is sensitive to spin and topol- 
ogy, the anomalous HaU effect (AHE)^ is at the focus of 
considerable theoretical and experimental efforts. The in- 
terest in spin sensitive phenomena is linked to the emerg- 
ing field of spintronics^, which offers an alternative to 
conventional charge-based electronics. The interest in 
the effects of topological features of bands on transport 
properties is linked to the role that these effects play in 
systems such as topological insulators and quantum Hall 
systems^. 

The AHE is described phenomenologically as trans- 
verse resistivity p^jf^^ or transverse conductivity cr^^^ 

linked to the intrinsic magnetization M of a conductor. 
Various models have been proposed: (a) The extrinsic 
model relates the AHE to antisymmetric scattering pro- 
cesses and it provides that 

p^f ^ = RsiioM^ (1) 

where Rs = apxx + bp^.^, ^'^d M± is the component of 
magnetization perpendicular to the film. The linear term 
in resistivity of Rs is attributed to skew scattering^ and 
it is expected to dominate in high conductivity regime 
{cxx > 10^ cm~^). The quadratic term is attributed 
to side jumps^ and it is expected to dominate in the good 
conductivity regime {axx ~ — 10^ fi"^ cm~^). (b) 
The intrinsic model known also as the Karplus-Luttinger 
model (K-L)^ or Berry phase model attributes the AHE 
to intrinsic topological properties of the band^,*^. Ac- 
cording to this model p^j^^ = p'^x'^xyiM) and it is ex- 
pected to dominate in the same regime as the side jump 
mechanism. In the poor conductivity regime {(Txx < 
lO'* ri^^ cm^^) a universal behavior axy ^ f^^^^'* has 
been observed experimentally^; however, a theoretical 
understanding is still lacking. Interestingly, a similar 
scaling is predicted for metals in the limit of strong scat- 
tering due to finite-lifetime disorder broadeningiS, and 



within a microscopic model accounting for fiuctuations 
of local orbital energiesii. 

SrRuOs has played a pivotal role in the study of the 
AHE and numerous attempts have been made to eluci- 
date its complicated behavior. Berry phase calculations 
which assume a temperature-dependent exchange gap 
that closes at Tc seemed to describe the data reasonably^. 
However, a test of this scenario that focused on the van- 
ishing point of the AHE found that it vanishes for a given 
film (whose resistivity and magnetization were varied by 
field) at a specific resistivity, and not at a specific mag- 
netization as one may expect from a scenario which at- 
tributes the vanishing signal of the AHE to the Berry 
phase contribution at a particular exchange splitting^^-. 
Mid infrared measurements suggest the applicability of 
the Berry phase scenario at energies above 200 meV 
while the dc limit is dominated by extrinsic scattering 
mechanisms^. 

By using SrRuOa films with a wide range of thicknesses 
that vary considerably in the temperature-dependence of 
their resistivity, we provide a compelling piece of evi- 
dence that resistivity, irrespective of its sources or nature 
(elastic or inelastic), determines the AHE of SrRuOa 
in the entire ferromagnetic phase. This observation 
strongly suggests that changes in Berry phase due to 
assumed temperature-dependent exchange splitting can- 
not explain the complicated temperature dependence of 
the AHE. We show that the side jumps mechanism com- 
bined with the Karplus-Luttinger (Berry phase) mech- 
anism that takes into account the scattering time may 
explain the observed behavior. 



II. SAMPLES AND EXPERIMENT 

Our samples are epitaxial thin films of SrRuOa grown 
on slightly miscut (^ 2°) substrates of SrTiOa by reac- 
tive electron beam evaporation. The films are untwinned 
orthorhombic single-crystals, with lattice parameters of 
a = 5.53A, b = 5.57A, and c = 7.85A. The films were 
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FIG. 1: (a) Remanent AHE resistivity (p^jf ) of 8 films vs 
temperature, (b) Tlie AHE coefficient (Ra) vs temperature 
derived from pxy^ and A4"x using Eq. f . Inset: scaling of 
the perpendicular magnetization (Afx) normalized by its low 
temperature value as a function of temperature for 5 films 
with thickness between 6 to 90 nm. 

patterned to allow transverse and longitudinal resistivity 
measurements, which were performed with a Quantum 
Design PPMS-9. The films exhibit exceptionally high re- 
sistivity ratio (up to 90) indicative of their high quality. 
The thinnest films (< 10 nm) exhibit lower resistivity 
ratio (> 5) which is still very high considering the en- 
hanced surface scattering. Magnetic characterization of 
the films was performed using a Quantum Design SQUID 
magnetometer (MPMS). 

Magnetic films may exhibit AHE if their magnetiza- 
tion has a component perpendicular to the film plane. As 
shape anisotropy favors in-plane magnetization, in many 
cases a perpendicular field should be applied in order to 
tilt the magnetization out of the plane. This may com- 
plicate the analysis since the applied field also induces 
ordinary Hah effect (OHE). 

SrRuOs films exhibit intrinsic uniaxial magnetocrys- 
talline anisotropy with an easy axis which varies with 
temperature between 45 degrees to the normal at Tc to 
30 degrees at 2 K— . Moreover, the remanent magnetiza- 
tion is stable and spontaneous breakdown into magnetic 
domains occurs only a few degrees below T^. These 
features enable direct measurement of zero field (rema- 
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FIG. 2: (a) Longitudinal resistivity {pxx) of the 8 films pre- 
sented in Fig. [T] vs temperature, (b) Rs from Fig. [TJb) vs 
resistivity {pxx)- 



nent) antisymmetric transverse resistivity which can be 
fully attributed to AHE, p^"^. 



III. EXPERIMENTAL RESULTS 

Fig. [Ha) shows temperature dependence of p'^y^ of 
8 different samples. To extract Rs based on Eq. 1 (see 
Fig. mb)), we divide p^y^ by Mj, (shown in the inset). 
We note that for the thickness range of our samples, M±_ 
is practically identical except for small deviations related 
to thickness-dependent Tc- This is expected as thickness- 
induced changes in magnetic properties were reported for 
films with thickness lower than 6 nrtii^. We note that 
while Rs of the various samples has general common fea- 
tures, the variations are considerable. In particular we 
note differences in the values of Rs at 2 K, in the location 
of the negative peak and in the temperature at which Rs 
changes its sign. The large spread in Rs seems to corre- 
late with changes in the resistivity of the films, strongly 
affected by film thickness (see Fig. Ella)). However, as 
seen in Fig. [UJb), the extracted Rs does not scale with 
Pxx- In particular, we note that the resistivity at which 
Rs changes its sign (po) varies between 105 /ifJ cm for a 
50 nm thick sample to 202 /ifi cm for a 6 nm thick sam- 
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FIG. 3: Rs normalized by its absolute maximum value (Rs*) 
vs pxx normalized by its value when Rs changes its sign (p*). 
Inset: Rb* vs (p*)^. 
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FIG. 4: Normalized resistivity (p*) vs temperature. Inset: 
film thickness (d) divided by its resistivity when Rs changes 
its sign (po) vs film thickness (d). 



pie. Does this observation exclude the scenario that Rs 
is determined by pxx in the entire ferromagnetic phase? 
- not necessarily. 

Fig. 13] shows that i?*, defined as Rs normahzed by 
its maximum absolute value, does scale with p*, defined 
as Pxx normalized by po- The scaling function has a 
quadratic dependence on p* in the low resistivity regime 
(see inset) and it reaches its negative peak for all samples 
at p* ^ 0.7. 

A possible explanation for the striking scaling is that 
Rs is determined by pxx and that it does vanish at the 
same intrinsic resistivity p™* for all samples, consistent 
with a previous report^^; however, there is a multiplica- 
tive factor 7 between the nominal resistivity and the in- 
trinsic resistivity, pxx = IP^xx- ^ trivial source for 7 is 
uncertainty in film thickness and in geometrical factors 
of the pattern. However, these sources alone cannot ac- 
count for the observed variations of order 2. Another 
potential source is dead layer o^^'^'' whose existence may 
affect considerably the calculated resistivity of ultrathin 
films. Assuming a dead layer of thickness 5, we would ex- 
pect J — d/{d — 6) and a linear dependence between d/pa 
and d, as observed in the inset of Fig. The linear fit is 
consistent with a dead layer scenario with 6^3 nm and 
Pq"* ^ 100 p,ft cm. The dead layer scenario also implies 
that the resistivity of the various samples at high tem- 
peratures is not different (as suggested by Fig. [SJa)) but 
quite similar (Fig. 2]). As the main difference between the 
films is in their thickness, the result supports the dead 
layer scenario as it is expected that at high temperatures, 
where the mean free path is small and bulk scattering is 
dominant, the resistivity of our films would be similar. 

The dead layer scenario implies the need to normalize 
Rs; However, its division by 7 does not scale the data 
along the y axis. Therefore, the normalization of Rs 
with its maximum absolute value merely indicates that 
for all films there is a single Rs{pxx) function up to a 



multiplicative factor. 

We note that the scaling is obtained for films that vary 
considerably in their thickness and residual resistivity; 
namely the same value of Rs is obtained for very differ- 
ent values of M±^ and for very different contributions to 
Pxx- Thus for instance, Rs attains its maximum value 
at T/Tc=0.47 for the 6 nm thick film and at T/Tc=0.63 
for the 90 nm thick film. At this temperature the mag- 
netization is 84 percent (77 percent) of its low tempera- 
ture value for the thin thick) film and the resistivity is 2 
times (32 times) larger than its low temperature value. 
Therefore, point defects, surface scattering, magnons and 
phonons have very different weights in the two cases. 

IV. THEORETICAL MODEL 

The low temperature dependence of Rg on is 
consistent with side jumps mechanism^ and with the 
Karplus-Luttingcr (K-L) or Berry phase mechanism^. 
However, whereas side jumps can explain the scaling with 
Pxx due to its insensitivity to the scattering potential, it 
cannot explain the non monotonic temperature depen- 
dence which includes a sign change at higher tempera- 
tures. On the other hand, attributing the non monotonic 
temperature dependence to K-L mechanism with temper- 
ature dependent exchange gap as suggested previously^, 
yields p^y^ = p'ix'^xy^^) with a comphcated depen- 
dence of Uxy on M which is inconsistent with the scaling 
which assumes linear dependence on M^. 

We now show that a combination of the side jumps 
mechanism and the K-L mechanism which considers the 
effect of scattering rate and its temperature dependence 
(without assuming any change in the band structure) is 
a possible scenario. A consideration of the scattering 
rate (1/t) effect on the transverse conductivity in the 
K-L mechanism is required in moderately good conduc- 
tors, where h/r is not negligible compared to the inter- 
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FIG. 5: AHE coefficient Rs as a function of resistivity pxx for 
a thick film (500 A). Tlie dasfied line is a fit to Eq. 8. 



band gap. The leading correction yields a decrease of 
axy as the resistivity increases, and thus a possible non- 
monotonic behavior of p^y 

Accounting for a finite r, the K-L contribution to the 
AHE resistivity from Kubo's formula^^ becomes 



P. 



K-L 



E<nk\vy \ 7nkXmk\vx \ nk>( f (Cn^k) — f (sr. 
n^m,k {i(err^,k-e,^^k)-\-^/r}{er^^k-err^,k} 



(2) 



where is the crystal volume, k is the quasi-momentum 
n, m are band indices associated with the eigenvalues of 
the perfect crystal Hamiltonian, Vx, Vy are the velocity 
operators, and /(e) is the Fermi-Dirac distribution. This 
contribution accounts only for the intrinsic part of the 
AHE, i.e. ignores the corrections to the scattering pro- 
cesses due to spin-orbit interaction. We consider a model 
in which the main contribution to the sum in Eq. [2]is due 
to two bands denoted as 1,2, where the Fermi level crosses 
the upper band while the lower band is fully occupied. 
We further assume that the dominant contribution arises 
from states with quasi-momentum in a set denoted as K, 
and that the energy gap for fc 6 K between non-occupied 
states in the upper level and occupied states in the lower 
level is approximately independent of quasi-momentum 
and takes the characteristic value of A. Under these as- 
sumptions, we obtain 



K—L _ 2 2f-( 
Pxy ~ ~Pxx^ ""l 



A 



A* 



{iA-h^/rjA {-iA-hS/r}A 



where A is defined as 



A = 



< lk\vJ2k >< 2k\vx\lk > 



3 ^ J^,^fUy 

Zn 



) (3) 



(4) 



which can be associated with a Berry's phase^. Since A 
is odd under time-reversal and hence purely imaginary, 
we get 



Pxy 



2 2 fc 

-Pxx^ n- 



2Im{A) 
.2 + (ft/r) 



(5) 



Considering spin orbit interaction (SOI), A is expected 
to be proportional to M±^; namely. 



Im{A) = aM± 



(6) 



where a is a constant. Previous reports indicate that the 
band structure in SrRuOa is temperature independenti^, 
thus, ferromagnetism in SrRuOs should be described in 
the local band model^S. Therefore considering Eqs. 5 
and 6 and preforming averaging yield: 



rxy 



where M± is the averaged magnetization in the sample. 
Finally, we note that within the same level of approxima- 
tion (i.e., leading order in the scattering potential), the 
side jumps contribution is additive to the K-L term^i. 
Thus, the AHE coefficient (R^) is given by a sum of the 
two contributions: 

— pIx ^2 ^ {H/tY ^^^^ ■ 

The first term is the K-L term (with all the constants and 
the minus sign included in _B), and the second term is the 
side jumps contribution. As _B, A and C are merely as- 
sociated with the band structure they are assumed to be 
constants. 1/r is assumed to be, as usual, proportional 
to pxx and the proportionality factor is estimated based 
on band calculations^^. Thus, the right hand side in Eq. 
8 is a function of pxx alone. 

Fig. [5] shows a fit of our data using Eq. 8 where the 
parameter C is limited to an interval which corresponds 
to a reasonable range of side jumps (0.1-10 A^). We 
obtain a good fit for side jumps in the range of 1-10 A, 
and A in the range of 0.07-0.2 eV. The value of A is in 
good agreement with the characteristic energy at which 
lTii[axy) has a peak, measured in the infrared regime for 
the low temperature limit^'^. The fit presented in Fig. [5] 
is for a side jump ~ 4 A and A ~ 0.13 eV. The fact that 
a similar temperature dependence of R^ is observed for 
other systems^'^ suggests that this scenario is relevant to 
other materials as well. 



V. CONCLUSIONS 

The scaling of Px^^ data with pxx in SRO films im- 
plies that the AHE coefficient is determined by the to- 
tal resistivity irrespective of the relative contributions of 
different scattering processes. To explain the scaling and 
the non monotonic behavior of the scaling function we 
present a scenario that attributes the observed behavior 
to two contributions: (a) side jumps mechanism and (b) 
K-L (Berry phase) mechanism including the effect of fi- 
nite scattering rates. In the limit of low resistivity, the 
two contributions have quadratic dependence on resistiv- 
ity with coefficients of opposite signs where that of K-L 
term is larger. As resistivity increases, the K-L term 
decays due to the effect of finite scattering rates which 
yields a sign change of Pxy^- 
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